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INTRODUCTION 

The discovery of superconductivity in alkali- metal-doped fullerenes^ KsCeo and 
RbsCgo has raised interesting questions about the electron- phonon coupling in such 
compounds and its interplay with Coulomb repulsion. Ceo is a highly symmetri- 
cal molecule i.e. it is a truncated icosahedron and its electronic lowest unoccupied 
molecular orbitals (LUMO) are threefold degenerate^'"^''^. They form a Ti^ repre- 
sentation of the icosahedral group Ih- Filling the LUMO in Cqq~ anions leads in a 
naive picture to narrow, partially filled bands in the bulk fuUerides. The bandwidth 
W is determined by the hopping between the Ceo molecules which are quite far apart 
and W ~ 0.5 eV. The coupling of some Hg phono ns with electrons residing in the 
Ti^ orbital has been suggested to be responsible for the superconductivity^'^'^. The 
Coulomb repulsion also may be important on the ball^. Several authors^"^^ have 
undertaken the study of the Jahn- Teller distortion that is expected in the fullerene 
anions. In such calculations one approach is to consider the electrons as fast degrees 
of freedom and the phonon normal coordinates are treated as static^^. This is the 
strong-coupling approach to the Jahn- Teller effect. Of course it is more difficult 
to treat the fully dynamical problem of the phonon mode coupled to the electrons. 
Here we will investigate the interplay between the electronic and phononic degrees 
of freedom on an isolated fullerene anion. We use a weak coupling approach to the 
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Jahn- Teller effect The lifting of degeneracy is obtained by a perturbation calcula- 
tion in the case of an undistorted anion. The ordering of levels can be described as 
" anti-Hund" ' rule. This calculation is very close in spirit to the standard treatment 
of the electron- phonon coupling in superconducting metals. In fact we show that 
this weak-coupling approach is exactly parallel to the nuclear physics calculation 
of pairing of nucleons in a single shell. We recover the so-called seniority model 
of pairing. It is often claimed that "superconductivity is a dynamical Jahn-TcUer 
effect" . While such a statement is not very appropriate to bulk metallic supercon- 
ductors, we will show that in the context of fuUerene anions, it really makes sense. 
There is now some evidence that the crossover from weak to strong coupling in 
this Jahn- Teller problem is smooth^^, so the weak-coupling calculation captures the 
physics of the intermediate regime. 

As in nuclear physics we find even-odd effects due to pairing in the ions. The 
effects we observe may be sought by spectroscopy of solutions of fuUerides in liquid 
ammonia, for example. We discuss the opposite effect of Coulomb interaction, lead- 
ing to Hund's rule in ordinary situations. Finally we point out that experimentally 
observed spectra may be at least partially explained by our calculation^^. 

THE ON-BALL ELECTRON-PHONON INTERACTION 

The electronic structure of tt electrons in the Cgo molecule is well known to 
be given by a simple Hiickel calculation. The levels are labeled^ by the irreducible 
representations (irreps) of the icosahedron group In- One important property has 
to be noted: three of the Ih irreps are the 1—0,1,2 spherical harmonics of SO (3) 
which do not split under the Ih group. They are commonly named A^, Ti„, H^. In 
addition there is also the twofold spin degeneracy. 

In the ground-state of the neutral Ceo molecule all levels up to H^j included 
are completely filled thus building a singlet state I^I/o)- The LUMO are the six Ti^j 
states. These are occupied upon doping with extra electrons and the ground-state 
becomes then degenerate. One then expects the Jahn-Teller effect to distort the 
anion and lift this orbital degeneracy^''' due to the coupling of the Ti„ electrons 
to the vibrational modes of the molecule (also referred to as phonons). In a weak- 
coupling scheme the phonon is purely virtual and is exchanged between electrons. 
Phonon exchange between electrons leads to an effective electron-electron interac- 
tion that competes with Coulomb repulsion and may lead to anti-Hund ordering of 
energy-levels. 

As a first investigation of electron-phonon coupling we use a perturbation 
scheme suited to degenerate levels we will derive an effective electron-electron inter- 
action with the assumption that filled states lying below the Ti^ level remain frozen 
so that intermediate states involve only Ti„-Ti^ excitations. Indeed the H^-Ti^( 
gap is fa 2eV whereas maximum phonon energies are 0.2eV. 

A typical electron-phonon interaction term reads: 




amim2 a miCT m2<T' 
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Here Xq, are normal coordinates, the subscript referring both to the irrep and to the 
row in the irrep they belong to, cj^ ^ is the creation operator for an electron with 
spin a in the Ti^ (^=1) level, mi taking one of the m=— 1,0,1 values, and famim2 
are complex coefficients. The c}^^ operators transform as 1=1 \l,m) vectors un- 
der Ih symmetries, and their conjugates Cma transform as —m) vectors. 
The (— l)"^^"'""^c]„io-C-m2cr products transform then as members of the Ti^ x Ti„ 
representation, which in the Ih group splits as: 

Ti„ X Ti„ = Ag + Tig + Hg. 

This selects the possible vibrational modes Ti„ electrons can couple to. In fact, 
only Hg modes split the degeneracy^. 

Let us consider a particular fivefold degenerate multiplet of modes. Their 
normal coordinates will be labelled X^, m ranging from -2 to +2. Since appears 
only once in the product Ti„ x Ti^ , the interaction is determined up to one coupling 
constant g by the usual formula for the coupling of two equal angular momenta to 
zero total angular momentum: 

W = gY,{-^TXm^-m. (1) 

m 

The may be chosen such that = {—1)'^X_^ and have the following expres- 
sion in terms of phonon operators: 

= ^ {am + i-^Tal^) (2) 

whereas the $rn are the irreducible 1=2 tensor operators built from the c^c products 
according to: 

= J](l, 1, 2|mi, m - mi, m)(-l)(— ^+^)4^,c_^+^^,, (3) 

mi 

where (ii, ^2, ^|"7,i, m2, m.) arc Clcbsch-Gordan coefficients. 

We now consider a doped CgQ~, molecule, < n < 6. Its unperturbed degener- 
ate ground-states consist of to which n Tiu electrons have been added times 
a zero- phonon state. They span a subspace denoted by £q. In £q the unperturbed 
Hamiltonian Hq reads: 

Ho = eti„ Yl ^^rnaCrna + ^1 ^m^m^ 
m,a m 

where et^^ is the energy of the Ti^t level, hu) is the phonon energy of the mul- 
tiplet under consideration. Within Sq the effective Hamiltonian up to second order 
perturbation theory is given by: 

Heff = EoPo + PoWPo + PoW{l - Po)—^-—{l - Po)WPo, 

^0 — -"0 
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where Pq is the projector onto Sq, Eq is the unperturbed energy in this subspace 
which is just the number of doping electrons times e^^^. The hnear term in W gives 
no contribution. Using expressions (1) and (2) for W and one finds: 

^ (-ir$™.,$-^.„ (4) 

m,CTlCT2 

where we have now included spin indices. We can now use equation (3) to express 
H(,ff as a function of c and c^^ operators and put it in normal ordered form using 
fermion anticommutation rules. In this process there appears a one-body interaction 
term which is a self-energy term. We will henceforth omit the Hq term which is a 
constant at fixed number of doping electrons. 

Let us now define pair creation operators ^f^^ which when operating on the 
vacuum |0) create pair states of Ti„ electrons that are eigenfunctions of L,S, Lz,Sz, 
where L,S are total angular momentum and spin, and Lz,Sz their ^-projections. I 
and s can take the values 0,1,2 and 0,1 respectively. This holds also if |0) is taken 
to be the singlet state l^'o)- 



mi,CTi 



The quantity Af^^ is non-zero only if (l + s) is even and the norm of ^f^^|0) is 
then equal to ^/2. The inverse formula expressing c^c^ products as A'^ operators is: 



cLiaiCLas = Xl^^' ^' "^1 + TO2)(-, -, s\ai, CT2, di + (^3) . (6) 

l,s 



As Heff is a scalar, its two-body part may be written as a linear combination of 
diagonal Af^'^Af^ products whose coefficients depend only on / and s: 

J2 ni.s)Az^Az. 

ls,ma 

The F{l,s) coefficients are calculated using expressions (4), (3), (6). We then get 
Hf,ff in final form: 



rr _ ^9 
^eff - 



6hu; 



I N -L. 400t .00 _ 1 4lo-t .la , J_ 4OO t .00 ) (Y) 

I -'V + ^00 ^00 2 1^ 10 ^ '^'^ I ' ^ ' 

\ m,a m / 



In this formula N is the electron number operator for the Ti„ level; the N term 
appears when bringing H^ff of expression (4) in normal ordered form. In our 
Hamiltonian formulation the effective interaction is instantaneous. 

There are actually eight multiplets in the vibrational spectrum of the Ceo 
molecule. To take all of them into account we only have to add up their respective 
coefficients 5g'^/6?iLU, their sum will be called A. 



4 



THE ELECTRONIC STATES OF FULLERENE ANIONS 

We shall now, for each value of n between 1 and 6, find the n-particle states and 
diagonalize He/f. The Hamiltonian to be diagonalized is that of equation (7) where 
the prefactor is replaced by —A. The invariance group of H^ff is Ih x SU(2). The 
n-particle states may be chosen to be eigenstatcs of L, S, L^, and we shall label 
the multiplets by {I, s) couples, in standard spectroscopic notation (^*+^L stands for 
(l,s)). The pair (Z, s) label SO(3) x SU(2) irreps which, as previously mentioned, 
remain irreducible under Ih x SU(2) as long as I doesn't exceed 2; for larger values 
of / SO(3) irreps split under Ih- Fortunately enough, the relevant values of I never 
exceed 2. Moreover given any value of n, (/, s) multiplets appear at most once so 
that the energies are straightforwardly found by taking the expectation value of the 
Hamiltonian in one of the multiplet states. The degeneracies of the levels will then 
be {21 + l)(2s + 1). We now proceed to the construction of the states. 

• n=l: There are six degenerate states cJ^^|\E'o) whose energy is —A. 

• n=2: There are 15 states, generated by applying Af^^ operators on l^'o)- 
There is one state, five states and nine states with energies -4A, -llA/5, 
-A. 

• n=3: There are 20 states. States of given l,m,s,a can be built by taking 
linear combinations of ^''"c't'l^'o) states according to: 

mi,(Ti 

These states belong to the following multiplets: ^P (E=-3A), (E= -9 A/5), 
(E=0). 

• n=4: There are 15 states, which are obtained by applying Af!^^ operators 
on ^gg^|*o). They are (E=-4A), (E= -llA/15), ^P (E=-A) 

• n=5: There are six ^P states which are 

cJ„^ylggUgg^|*o) and whose energy 

is -A. 

• n=6: There is one state whose energy is 0. 

It is interesting to note that the above treatment of electron-phonon interaction 
parallels that of pairing forces in atomic nuclei^^'^^. Of course in the case of finite 
fermionic systems there is no breakdown of electron number but there arc well- 
known "odd-even" effects that appear in the spectrum. In our case pairing shows 
up in the ground state for CgQ rather than ^P as would be preferred by Coulomb 
repulsion i.e. Hund's rule. The construction of the states above is that of the 
seniority scheme in nuclear physics^^. We note that similar ideas have been put 
forward by V. Kresin some time ago, also in a molecular context^^. The effective 
interaction that he considered was induced by a core polarization. 
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THE EFFECT OF COULOMB REPULSION 



We now consider the Coulomb electron-electron interaction and assume it to 
be small enough so that it may be treated in perturbation theory. To get some 
feeling of the order of magnitude of this repulsion we use the limiting case of on-site 
interaction i.e. the Hubbard model. This Hamiltonian is not specially realistic but 
should contain some of the Hund's rule physics. The two-body interaction now 
reads: 

Usr^ t t 

2 / J ^ia^i—a^i—a^icn 
i,a 

where the i subscript now labels the tt orbitals on the Ceo molecule. The quantity U 
is ~ 2-3 eV from quantum chemistry calculations^^ Since level degeneracies are split 
at first order in perturbation theory we confine our calculation to this order and 
have thus to diagonalize the perturbation within the same subspace So as before. 
In this subspace it reads: 

Wh = U {a\i){P\i)mm ci^cl^c^^csp 

i, a/375 

where greek indices label one-particle states belonging either to \^o) or to the T^y^ 
level. Let us review the different parts of Wh- Note that since the j^'o) singlet 
remains frozen we have the identity: cjj,c^ = 5q;/3 if ct, P label states belonging to 
l*o). 

-A part involving states belonging to j^^o) only: 

i,al3 

a,(3 belong to l^'o)- This term is thus diagonal within £q and merely shifts the total 
energy by a constant that does not depend on the number of doping electrons. It 
won't be considered in the following. 

-A part involving both states belonging to l^'o) and to the Ti„ level: 
where a, 6 belong to the Ti^ level whereas /? belongs to l^'o)- It reduces to: 

aS,a ^ i (3 ^ 

The sum over /3 is just the density on site i for a given spin direction of all states 
belonging to |\E'o) which is built out of completely filled irreps. As a result this 
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density is uniform and since j^^o) contains 30 electrons for each spin direction it is 
equal to 1/2. then becomes diagonal and reads: 

U x-^ t 

Its contribution is thus proportional to the number of Ti^j electrons. It represents 
the interaction of the latter with those of the singlet and we won't consider it in 
the following. 

-A part involving only states belonging to the Ti„ level: 
has the same form as Wh with all indices now belonging to the Ti„ level. 
Whereas the interaction has a simple expression in the basis of \i) states, we need 
its matrix elements in the basis of the Ti„ states. There are in fact two Ti^ triplets 
in the one-particle spectrum of the Cgo molecule, the one under consideration having 
higher energy. To construct the latter we have first constructed two independent 
sets of states which transform as x,y, z under Ih. These are given by: 



ex. 



= E 



and 



\a 



where Cq, are three orthonormal vectors, i labels sites on the molecule, the fi are the 
vectors joining the center of the molecule to the sites while the ki join the centre of 
the pentagonal face of the molecule the site i belongs to to the site i. We assume 
that the bonds all have the same length. These states span the space of the two Ti^ 
triplets. The diagonalization of the tight-binding Hamiltonian in the subspace of 
these six vectors yields then the right linear combination of the |a) and \a)' states 
for the upper lying triplet. From the x, y, z states one constructs 1=1 spherical 
harmonics. We then get the matrix elements of Whs in the basis of Ti„ states. As 
So is invariant under Ih operations and spin rotations, Wh^ which is the restriction 
of Wh to So is invariant too. It may thus be expressed using the A, A'' operators by 
using formula (5) in the same way as the phonon-driven interaction and we finally 
get: 

Whs = (^^A'o'a'o'o + 3^ E ^Z'aZ^ , (8) 

which is the only part in Wh that we will keep. Note that there is no contribution 
from /=1, s=l A'^A products. Indeed the Hubbard interaction is invariant under 
spin rotation and couples electrons having zero total S^- As the coefficients of A^ A 
products depend solely on I and s they must be zero for s 0. The spectrum 
for any number of Ti^ electrons is now easily found. Of course the order of the 
multiplet is now reversed: for n=2, we have "^P, then ^D, then^S. For n=3, we have 
^S, then ^D, then ^P. For n=4, we have ^P, then ^D, then ^S. 
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CONCLUSION 



The ordering of energy levels in the electron-phonon scheme are clearly opposite 
to those of Hund's rule. The clear signature of what we can call "on-ball" pairing 
is the ground state of CqqI the two extra electrons are paired by the electron- 
phonon coupling. We note that the U of the Hubbard model appears divided by 
large factors: this is simply due to the fact that the Ceo molecule is large. As 
a consequence, if U ~ 2 eV, Coulomb repulsion may be overwhelmed by phonon 
exchange. With a H^, phonon of typical energy 100 meV and coupling 0(1) as 
suggested by numerous calculations^'^'^'^, the quantity A may be tens of meV. 

It seems to us that the cleanest way to probe this intramolecular pairing would 
be to look at solutions of fuUerides leading to free anions such as liquid ammo- 
nia solutions or organic solvents"*^^"^^. EPR or IR spectroscopy should be able to 
discriminate between the two types of spectra. Measurements by EPR should de- 
termine whether or not the two extra electrons in Cqq are paired, for example. 
In near-IR spectroscopy the lowest allowed transition for Cqq should be at higher 
energy than that of Cgg due to the pairing energy while in the Coulomb-Hubbard 
case it is at lower energy. 

Present experiments^^'^° have studied the near-IR spectra of solutions of ful- 
leride anions prepared by electrochemical reduction. There are several peaks that 
do not fit a simple Hiickel scheme of levels. They do not have an immediate inter- 
pretation in terms of vibrational structure^^'^'^. With our energy levels in table I, 
a tentative fit would lead to A f« 80 meV assuming U = 0. Such a value leads to 
intriguing agreement with the major peaks seen for Cgg and CgQ while this is no 
longer the case for CgQ and CgQ . 

Finally we mention that recent EPR experiments^^ have given some evidence 
for non-Hund behaviour of the fuUeride anions. While one may observe some trends 
similar to the results of the phonon-exchange approximation, it is clear that the 
model we used is very crude. In a bulk conducting solid we do not expect the 
previous scheme to be valid since the levels are broadened into bands: then phonon 
exchange leads of course to superconductivity. 
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